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Abstract 

The effects of interfacial roughness and oxide film thickness on thermally-induced stresses 
in plasma-sprayed thermal barrier coatings subjected to thermal cycling are investigated using 
the recently developed higher-order theory for functionally graded materials. The higher-order 
theory is shown to be a viable alternative to the finite-element approach, capable of modeling 
different interfacial roughness architectures in the presence of an aluminum oxide layer and 
capturing the high stress gradients that occur at the top coat/bond coat interface. The oxide 
layer thickness is demonstrated to have a substantially greater effect on the evolution of residual 
stresses than local variations in interfacial roughness. Further, the location of delamination 
initiation in the top coat is predicted to change with increasing oxide layer thickness. This result 
can be used to optimize the thickness of a pre-oxidized layer introduced at the top coat/bond coat 
interface in order to enhance TBC durability as suggested by some researchers. The results of our 
investigation also support a recently proposed hypothesis regarding delamination initiation and 
propagation in the presence of an evolving bond coat oxidation, while pointing to the importance 
of interfacial roughness details and specimen geometry in modeling this phenomenon. 


1 INTRODUCTION 

Many factors may affect the durability of ceramic thermal barrier coatings (TBCs). These 
include the thermoelastic and viscoplastic parameters controlling the response of the individual 
constituents (i.e., the top coat, bond coat and substrate, see Fig. 1), bond coat/top coat interface 
roughness and strength, bond coat oxidation, top coat microstructure (which depends on the fab- 
rication process [1]), substrate curvature, and imposed thermal load profiles. Each of these factors 
influences in a different manner the internal stresses induced in the individual constituents during 
service conditions, leading to potentially different failure modes [2,3]. A thorough understanding of 
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the influence of each of these factors separately, as well as their interaction, is therefore essential in 
developing guidelines for the design of durable TBCs. This requires the development of sufficiently 
comprehensive analytical and/or numerical models capable of accurately assessing the importance 
of the above factors on the evolution of internal stresses during service conditions. 

Early experimental investigations into failure mechanisms of plasma-sprayed TBCs conducted 
by Miller and Lowell [4] identified delamination at (or just above) the top coat /bond coat interface 
as the primary failure mode, ultimately leading to spallation of the coating. The failure mecha- 
nism was hypothesized to involve residual stresses induced into the coating during the cooldown 
stage of a thermal cycle due to the thermal expansion mismatch between the top coat and the 
bond coat. The initiation of the failure process was shown to be substantially influenced by the 
formation of a thin oxide film at the top coat/bond coat interface during exposure to elevated 
temperatures, and based on the observed location of delamination the interfaced irregularity was 
thought to play a role. In order to better understand the influence of bond coat/top coat rough- 
ness and bond coat oxidation on thermally-induced stresses in plasma-sprayed TBCs, and thus the 
observed spallation-type failure modes, finite-element investigations were conducted by Chang et 
al. [5, 6, 7, 8] on a model cylindrical TBC geometry with a top coat/bond coat interface varying 
sinusoidally in the circumferential direction at a fixed radius from the cylinder’s center. These 
studies demonstrated the evolution of a significant radial stress component at the top coat/bond 
coat interface upon cooldown, which became tensile at the wavy interface’s crest. Incorporation 
of plasticity into the finite-element analysis decreased the magnitude of the radial stress, while 
introduction of an oxide layer between the top coat and the bond coat produced a change in the 
sign of this stress along the oxide layer/top coat interface. Based on the finite-element results, a 
model was proposed for oxidation-induced coating failure in plasma-sprayed coatings with rough 
interfaces. However, validation of the proposed model required the incorporation of all the features 
investigated separately into a unified model, including more realistic thermal loading, appropriate 
inelastic models for the constituent phases, oxide layer evolution and crack growth. 

The results of similar, but more realistic, finite-element analyses, addressing the influence of 
bond coat/top coat interface roughness and bond coat oxidation on thermally induced stresses in 
plasma-sprayed ceramic coatings on cylindrical specimens, were recently reported in a sequence of 
papers by Ferguson et al. [9], Petrus and Ferguson [10] and FVeborg et al. [11]. In these studies, 
cyclic thermal loading that closely simulated actual in-service thermal exposure was employed, 
oxide layer growth was modeled explicitly, and a creep model was used for the inelastic behavior of 
the TBC constituents to account for the actual time-dependent response at elevated temperatures. 
The results indicate the evolution of a non-zero stress component normal to the interface during 
cyclic thermal loading. Due to stress relaxation during the hold period of a heatup-hold-cooldown 
thermal cycle, this stress component becomes tensile at the interface’s crest upon cooldown. The 
presence of this stress component correlates well with the observed horizontal cracks at or just above 
the interfacial crests which have been demonstrated previously to produce spallation-type failures 
in plasma-sprayed TBCs. Oxidation of the bond coat was shown to affect the stress distribution 
normal to the bond coat /top coat interface in a manner that accelerates the spallation process. 

In this paper, we report the results of an investigation which addresses similar issues but from 
a somewhat different perspective. First, we consider a flat TBC system shown in Fig. 1(a) that 
remains in a generalized plane strain state under the imposed spatially uniform cyclic thermal load- 
ing which simulates uniformly-applied furnace thermal cycling. Experiments on a similar geometry 
under this type of thermal loading were recently conducted by Koolloos and Houben [12]. Thus we 
eliminate any effects due to curvature that may cause additional radial stresses normal to the bond 
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coat/top coat interface. This allows us to investigate the stresses induced strictly by the interfacial 
waviness. Second, we consider a periodically varying wavy interface between the top and bond 
coats w'hose waveform consists of two sine waves which are mirror images of each other as shown 
in Fig. 1(a), in contrast with the single sine wave representation employed in the finite-element 
studies discussed above (i.e., Ref. [5-11]). The double sine wave contains cusps and thus is a more 
realistic representation of the locally irregular interfacial features produced by the plasma-spray 
process than the smooth single sine wave representation [3]. To investigate the effect of interfacial 
roughness on the stress field during the applied thermal cycle, the trough amplitude is varied (with 
the crest amplitude fixed) to produce low, standard and high degrees of interfacial roughness. The 
chosen model geometries allow one to investigate disturbances in the stress field produced by local 
interfacial roughness irregularities observed in plasma-sprayed coatings. Due to the periodic nature 
of the TBC system’s architecture, the analysis of thermally-induced stress fields is performed on the 
repeating cross-section element shown in Fig. 1 (b) under constraints which simulate generalized 
plane strain deformation and take advantage of the symmetry considerations. 

In addition, an oxide layer is introduced between the top coat and the bond coat of the standard 
interface roughness geometry to investigate its effect on the stress field. In this investigation, we 
do not attempt to simulate the growth of this oxide layer as a function of time as will be discussed 
in more detail later, in contrast with the work of Ferguson et al. [9], Petrus and Ferguson [10], and 
Freborg et al. [11]. Rather, we investigate how oxide layers with different thicknesses affect the 
thermally-induced stress field, with the mean thickness of the bond coat kept fixed. This simulates 
the introduction of a pre-oxidized layer between the top and bond coats to retard or control the 
rate of oxygen diffusion into the metallic bond coat, as was done by Koolloos and Houben [12], and 
allows one to quantify the oxide layer’s thickness effect on the stress field. 

The analysis of the thermally-induced stresses in the above configurations is conducted using 
the recently developed version of the higher-order theory for bidirectionally functionally graded 
materials, referred to as HOTFGM-2D in the literature [13,14,15,16]. As described in the next 
section, this analysis technique offers a comprehensive approach towards modeling the response of 
material systems with different microstructural details, including certain advantages not available 
in standard displacement-based finite element techniques. In particular, it is capable of capturing 
different microstructural details observed in plasma-sprayed coatings, including top coat’s porosity 
and oriented microstructure whose effects remain to be addressed in the literature. The predictive 
capability of this novel analytical technique demonstrated in this investigation sets the stage for 
future investigations that will address the effect of these (as well as other) microstructural details 
on TBC failure mechanisms and durability. 

2 MODELING APPROACH 

HOTFGM-2D is based on a geometric model of a heterogeneous composite having finite 
dimensions in the x 2 — 23 plane and extending to infinity in the x\ direction, Fig. 2. In the 22 — £3 
plane, the composite is functionally graded with an arbitrary distribution of inclusions of arbitrary 
cross section. The loading applied to the composite may involve an arbitrary temperature or heat 
flux distribution and mechanical effects represented by a combination of surface displacements 
and/or tractions in the 22 - 23 plane, and a uniform strain in the x\ direction. 

The functionally graded microstructure in the 22 — 23 plane is modeled by discretizing the 
heterogeneous composite’s cross section into N q and N T generic cells in the intervals 0 < 22 < if, 
0 < 23 < L. The generic cell (q, r) used to construct the composite, highlighted in Fig. 2, consists 
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of four subcells designated by the pair ($7), where each index /?, 7 takes on the values 1 or 2 which 
indicate the relative position of the given subcell along the £2 and £3 axis, respectively. The indices 
q and r, whose ranges are q = 1, 2, N g and r = 1,2, N r , identify the generic cell in the x 2 - £3 
plane. The dimensions of the generic cell along the x 2 and £3 axes, hff\ h^\ and l%\ can vary 
in an arbitrary fashion such that 


H = Y,(h? + h^) I = £((M + !«) (D 

<?=1 r=l 


Given the applied thermomechanical loading, an approximate solution for the temperature and 
displacement fields is constructed based on volumetric averaging of the field equations together with 
the imposition of boundary and continuity conditions in an average sense between the subvolumes 
used to characterize the material’s microstructure. As described briefly in the appendix, this is 
accomplished by approximating the temperature and inplane displacement fields in each subcell 
of a generic cell using a quadratic expansion in the local coordinates xy, x^ centered at the 
subcell’s centroid. A higher-order representation of the temperature and displacement fields is 
necessary in order to capture the local effects created by the thermomechanical field gradients, the 
microstructure of the composite and the finite dimensions in the functionally graded directions. 
The out-of-plane displacement, on the other hand, is approximated using linear expansion in local 
coordinates which reflects the material’s continuous microstructure in the out-of-plane direction. 
The unknown coefficients associated with each term in the temperature and displacement field 
expansions are then obtained by constructing systems of equations that satisfy the requirements 
of a standard boundary-value problem for the given field variable approximations. That is, the 
heat and equilibrium equations are satisfied in a volumetric sense, and the thermal and heat flux 
continuity conditions, and the displacement and traction continuity conditions, within a given cell, 
as well as between a given cell and its adjacent neighbors, are imposed in an average sense. 

Similarities and differences between this approach and standard displacement-based finite- 
element formulations have been described elsewhere [17]. An important advantage of the present 
approach is the simultaneous satisfaction of displacement and traction continuity conditions be- 
tween the different subvolumes of the spatially variable microstructure. This facilitates the de- 
termination of thermally-induced stresses across the top coat/bond coat interface in the present 
application which cannot be easily captured due to very localized and very steep gradients in the 
thickness direction as will be demonstrated. 

The volume discretization employed in HOTFGM-2D allows one to represent the geometries of 
different types of spatially variable material architectures in sufficient detail. For instance, this ap- 
proach was recently employed in modeling thermally-induced stresses in functionally graded TBCs 
subjected to thermal gradient cyclic loading, illustrating how creep-induced stress redistribution 
depends on the level of microstructural refinement [17]. In the present application of HOTFGM- 
2D, we discretize the repeating cross-section element of the TBC system shown in Fig. 1(b) in 
order to model the geometric details associated with different degrees of interfacial roughness and 
different oxide layer thickness. The mean top coat and bond coat thickness are the same as those 
investigated by Freborg et al. [11], using circular specimen geometry, while the substrate thickness 
is different due to the absence of curvature in our specimen geometry. The results of this volume 
discretization in the region around the wavy interface are shown in Fig. 3 for the three levels of 
interfacial roughness, herein called low, standard and high. The wavelength of a single sinusoidal 
waveform is 50 fim and the crest amplitude is fixed at 10 /im while the trough amplitude is 4, 10 
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and 16 fim for the low, standard and high interfacial roughness, respectively. The standard inter- 
face roughness in our investigation has the same wavelength to amplitude ratio as that employed 
by Ferguson et al. [9]. Included in Fig. 3 are the dimensions of thin and thick oxide layers for the 
standard interfacial roughness. As is observed, the chosen volume discretization adequately depicts 
the different levels of interfacial roughness. 

The material thermoelastic and viscoplastic parameters that govern the response of the indi- 
vidual TBC constituents are given in Table 1 [9]. The viscoplastic response of the Zr-based top 
coat, NiCr-based bond coat, and Ni-based substrate is modeled by a power-law creep equation, 
generalized to multi-axial loading situations as follows, 

_ 3F(q e ,r) , 

« 2u e ij 


where d' l3 are the components of the stress deviator, a e = yj3/2a^a^, and 


F(a e ,T) = A(T)a^ 


(3) 


where the temperature-dependent parameters A, n are listed in the table in the range 10 — 1200° C. 
These are the same as the parameters used by Freborg et al. [11], as are the temperature-dependent 
thermal expansion coefficients. However, average values of the Young’s moduli over the specified 
temperature range provided by the above authors were used in our calculations for computational 
efficiency as discussed below. 


Table 1. Thermoelastic and creep material parameters of the TBC constituents 

(after Ferguson et al. [9]). 


Material 

E (GPa) 

V 

a ( 10 " B /°C ) 
10 ° -> 1200 ° C 

A (MPa -n /s) 
10° — 1200 °C 

n 

10 ° -> 1200 ° C 

Zr-based TC 

25.4 

0.33 

7.6 -v 12.5 

2.01 x 10 _a ° -» 1.85 x 10 " 7 

1.0 -» 3.0 

AI 3 O 2 layer 

380 

0.26 

8.6 

— 

- 

NiCr-based BC 

156 

0.27 

12.0 -» 14.4 

4.39 x 10~ 4 ° -*■ 7.40 x 10 " 6 

O 

CO 

T 

O 

i-H 

Ni-based SUB 

156 

0.27 

12.0 19.3 

4.85 x 10 ~ 36 -*• 2.25 x 10 ~ 9 

1.0 -» 3.0 


The creep model employed in this investigation requires considerable computational effort to 
solve the governing systems of equations described in the appendix which characterize the tem- 
perature and deformation fields in the individual subcells of the discretized repeating cross section 
element. This is due to the high creep, and thus relaxation, rates at elevated temperatures that the 
TBC constituents experience for the given creep parameters and thermal loading. In the present 
investigation, a forward Euler integration technique was employed to integrate the governing equa- 
tions using a variable time step. Very small time steps were required during the heatup and 
cooldown portions of the thermal cycle to obtain convergent solutions. This, in turn, necessitated 
the use of average Young’s modulus values over the considered temperature range in order to avoid 
the need to invert the large global stiffness matrix in the governing equations at each time incre- 
ment. Convergent solutions were obtained as verified by repeating the calculations for decreasingly 
smaller time increments, as well as by ensuring that the continuity of interfacial tractions across 
the top coat /bond coat interface was satisfied. 
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Finally, the spatially-uniform thermal loading applied to the TBC system is shown in Fig. 4. 
It consists of a heatup, hold, and cooldown portions over a period of 630 sec. Again, this is the 
same loading cycle as that employed by Freborg et al [11]. Three thermal cycles were sufficient 
to delineate the major response characteristics for the considered material combinations, geometry 
and temperature range due to the rapid stress relaxation within the TBC constituents during 
the hold period. We note that transient effects were not taken into account in performing the 
calculations (i.e., steady-state thermomechanical analysis was employed) due to the relatively slow 
rate of heating/cooling employed in the thermal loading history (which, as previously mentioned, 
simulates furnace heating). Transient effects would certainly become important when simulating 
thermal shock tests, and can easily be accommodated by the higher-order theory. In addition, 
residual stresses were not included in the calculations due to the complete stress relaxation that 
occurs during the hold portion of the thermal cycle for the employed thermal history, as will be 
seen in the next section. It must be mentioned, however, that for other thermal histories for 
which complete stress relaxation does not occur in the individual constituents (in the presence of 
a through-thickness temperature gradient, for example), fabrication-induced residual stresses may 
play a role. Further, residual stresses due to fabrication may induce initial damage which may 
accelerate the spallation mechanism. It has been reported, for instance, that spallation occurs 
sooner when a coating is plasma-sprayed onto an excessively hot substrate [ 18 ]. 

3 RESULTS 

The results presented herein emphasize the influence of interfacial roughness and fixed oxide 
layer thickness on the through-thickness (<733) stress that is thought to be responsible for the 
formation of horizontal cracks in plasma-sprayed coatings in the vicinity of the top coat/bond coat 
interface’s crest. However, the transverse (1J22) stress is also included in the section addressing the 
oxide layer thickness effect as large magnitudes of this stress arise in the oxide layer itself and may 
potentially lead to the oxide film’s degradation. These stresses are presented at critical locations as 
a function of time, as well as through critical cross sections (both vertical and horizontal) at those 
points in time that correspond to the maximum magnitudes of these stress components. 

The presentation of the results is divided into two parts. In the first part, we initially demon- 
strate the convergence of our solution technique by focusing on the interfacial traction continuity 
across the top coat/bond coat interface, and then investigate the influence of interfacial roughness 
on the through-thickness stress in the absence of an oxide layer. In the second part, we investigate 
the effect of varying the oxide layer thickness on both stress components for the standard interfacial 
roughness, and compare the results with those when the oxide layer is absent. 

3.1 Interface Roughness Effects 

Figure 5 illustrates the evolution of the through-thickness stress across the crest and trough 
of the top coat/bond coat interface during the first thermal cycle for a low and a high interface 
roughness in the absence of an oxide layer. The stresses in the center of two adjacent subcells lying 
on the top coat and bond coat side of the interface’s crest and trough are displayed in order to 
demonstrate the interfacial traction continuity satisfied by the higher-order theory formulation. As 
mentioned in the foregoing, the satisfaction of the interfacial traction continuity in the presence 
of very rapid relaxation of the stress field at elevated temperatures requires careful integration 
of the governing systems of equations presented in the appendix. The demonstrated traction 
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continuity gives us confidence in the results that will be discussed in the sequel. To achieve a 
comparable level of traction continuity satisfaction using the displacement-based finite-element 
approach requires substantial mesh refinement in the interfacial region where the stress gradients 
become very pronounced during heatup and cooldown. The level of traction continuity satisfaction 
observed in our results may be difficult to achieve using the displacement-based finite-element 
analysis as suggested by the initial results reported by Ferguson et al. [9]. 

The stress profiles in Fig. 5 indicate that the maximum magnitudes of the through-thickness 
stress occur during the initial heatup portion of the thermal cycle (at about 30 sec) and after 
cooldown (at 630 sec). During the initial heatup period, the through-thickness stress becomes 
compressive across the interface’s crest and tensile across the interface’s trough. After the stresses 
at these locations reach their maximum values at 30 sec or 825° C, they quickly decay to zero 
before the initial heatup is completed due to the rapid stress relaxation at the elevated temperature. 
Therefore, during the cooldown portion of the thermal cycle, the through-thickness stress becomes 
tensile across the interface’s crest and compressive across the interface’s trough, with the maximum 
magnitudes occurring upon the cycle’s completion at 630 sec. During subsequent thermal cycles 
(not showm), the through-thickness stress across the interfacial crest and trough initially decreases 
to zero during the heatup portion and subsequently retraces its path during cooldowm, attaining 
the same maximum values as those after the first cycle. Thus the trough, in the absence of an oxide 
layer, will be subjected to either zero or compressive through-thickness stress during subsequent 
thermal cycling for the chosen TBC material system and thermal loading profile, eliminating the 
possibility of horizontal crack growth in this region. This will be true whether or not the interface at 
the trough survives the initial heatup. Alternatively, the crest will be subjected to tensile through- 
thickness stress during subsequent thermal cycles upon cooldowm, which provides the driving force 
for horizontal crack growth initiated by local delamination. 

The results of Fig. 5 indicate that varying the interfacial roughness by keeping the interfacial 
crest amplitude fixed while changing the interfacial trough amplitude has a substantial influence 
on the through-thickness stress across the trough, wnth the stress across the crest affected to a very 
small extent. This localized effect is more fully illustrated in Fig. 6 which shows the through- 
thickness stress evolution across the interface’s crest and trough for the three levels of interfacial 
roughness. At the crest location, the maximum through-thickness stress upon completion of the 
first and subsequent cycles occurs at the bond coat/top coat interface for all interface roughnesses. 
For the standard interface roughness it is 65 MPa. Decreasing the interface roughness reduces the 
maximum stress to 59 MPa, and increasing the interface roughness increases the maximum stress 
to 69 MPa. At the trough location, the maximum through-thickness stress occurs at the bond 
coat /top coat interface only for the high interface roughness as seen in Fig. 7 (see Fig. 3(a) for 
interface locations), which shows the stress distributions along the vertical cross section passing 
through the interface’s trough at the end of the first and subsequent cycles. The maximum through- 
thickness stress for the high interface roughness is —153 MPa. For the standard interface roughness 
the maximum stress is —99 MPa and occurs in the bond coat 3 [im below the bond coat/top coat 
interface, decreasing to —80 MPa at the interface. The maximum through-thickness stress for the 
low interface roughness, —55 MPa, also occurs in the bond coat, 3 \im below the bond coat/top 
coat interface, decreasing to —47 MPa directly at the interface. 

The through-thickness stress magnitudes at the crest and trough of the bond coat/top coat 
interface of the standard interface roughness generated by our analysis are comparable to the 
results obtained by Freborg et al. [11]. These authors have obtained 75 and -50 MPa at the crest 
and trough locations upon cooldown, respectively, versus 65 and —80 MPa obtained by our analysis. 
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The observed differences are likely due to the difference in the interfacial roughness representation 
(i.e. , single vs double sine waveforms) and specimen geometry (curved vs flat geometry). 

It is important to note that if the manner of varying the interface roughness were reversed 
by keeping the trough amplitude fixed while changing the crest amplitude, a smilar effect will 
also be observed. That is, the interfacial through-thickness stress at the trough would be affected 
very little, with the stress at the crest affected to a similar extent as at the trough in the present 
situation. 

We end this section by presenting the horizontal (along the direction) through-thickness stress 
distributions at the elevation passing through the crest of the bond coat /top coat interface (23 = 640 
fj-m, see Fig. 3(a)) in Fig. 8 to demonstrate how they are affected by the interface roughness, and 
to provide comparison with the results of the following section. As observed (and suggested by the 
vertical distributions in Fig. 7), the tensile through-thickness stress at the crest location changes 
to compression in the top coat region above the trough. The low interface roughness decreases 
the magnitude of the compressive stress somewhat while the high interface roughness increases it 
relative to the stress distribution for the standard interface roughness. Corresponding distributions 
will be presented in the following section to demonstrate how the through-thickness stress to either 
side of the crest is dramatically affected by the introduction of an oxide layer. 

3.2 Oxide Layer Thickness Effects 

Next, we demonstrate the effect of oxide layer thickness on the through-thickness stress for 
the standard interface roughness, and then focus on the transverse stress. Figure 9 compares the 
through-thickness stress evolution with time at points just below the crest and trough locations (on 
the bond coat side of the bond coat/top coat and bond coat/oxide layer interface) of the interfaces 
without and with an oxide layer. At this particular location, the presence of the oxide layer increases 
the magnitude of the through-thickness stress during the intial heatup and subsequent cooldown to 
a substantially greater degree than the corresponding increase due to the interfacial roughness. As 
before, during subsequent cycles, the residual stresses decrease to zero during the heatup period, 
remain zero during the hold period, and retrace the same path during cooldown as in the first cycle. 

More insight into the oxide layer-induced stress amplification is obtained by considering the 
through-thickness stress distributions in the vertical cross sections passing through the crest and 
trough at the end of each cycle. These are presented in Fig. 10. In the cross section passing 
through the crest, the maximum stress (65 MPa) in the absence of the oxide layer occurs exactly 
at the bond coat/top coat interface (23 = 640 fim) as discussed previously. The introduction of a 
thin oxide layer increases the maximum stress to 169 MPA and moves its location 1 nm below the 
bond coat/oxide layer interface. At the bond coat/oxide layer and oxide layer/top coat interface 
the through-thickness stress rapidly decreases to 71 and 42 MPa, respectively, increasing again to 
52 MPa 2 fim above the oxide layer/top coat interface (#3 = 643 fxm), and subsequently decaying 
to zero in a manner similar to that when the oxide layer is absent. In the presence of the thick 
oxide layer the maximum stress further increases to 500 MPa and its location moves back up to the 
bond coat/oxide layer interface. The through-thickness stress rapidly decreases within the oxide 
layer with increasing distance from the bond coat/oxide layer interface, changing from tension to 
compression within the oxide layer itself in contrast to the thin layer case. It attains its greatest 
compressive value of 52 MPa at the oxide layer/top coat interface ( 3:3 = 643 fim) at which point it 
reverses direction, becoming tensile a small distance from this interface (21 MPa at = 650 ^m), 
reversing direction again and finally decaying to zero. 
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Similar qualitative behavior is observed in the cross section passing through the trough, except 
for the sign of the through-thickness stress. In this case, the location of the maximum magnitude 
of the compressive through-thickness stress relative to the bond coat interface is slightly different 
in each case. In the absence of the oxide layer, the compressive through-thickness stress attains 
its greatest magnitude at an elevation 3 fim below the bond coat/top coat interface, as discussed 
previously. The introduction of a thin oxide layer moves the maximum stress location to the bond 
coat/oxide layer interface. This is further moved up into the oxide layer itself 1 [xvn above the bond 
coat /oxide layer interface when the oxide layer thickness is increased. The decay of the through- 
thickness stress with increasing distance from the location where its magnitude is maximum is 
monotonic when the layer is thin, and has the same character as in the cross section passing 
through the crest (except for the sign). In the case of the thick oxide layer, however, the through- 
thickness stress reverses sign, becoming tensile outside the oxide layer, and attains a maximum 
value of 16 MPa at the elevation of 630 \xm (average bond coat/top coat interface) above which it 
decreases to zero. This is similar to the sign reversal in the through-thickness stress observed in the 
cross section passing through the crest for the thick oxide layer case, except now the sign reversal 
occurs in the top coat above the oxide layer/top coat interface. 

The substantial stress amplification due to the oxide layer suggests that the initiation of delam- 
ination at the crest upon cooldown may potentially occur at the bond coat/oxide layer interface 
for both oxide layer thicknesses if the interfacial bond is weak (in the case of the thin oxide layer, 
the delamination initiation in the bond coat just below the bond coat/oxide layer interface where 
the maximum stress occurs can probably be ruled out for a dense Ni-based alloy). Subsequent 
delamination growth along this crest elevation in the horizontal cross section depends on several 
factors, including the strength of the bond coat/oxide layer interface, the bond coat/oxide layer 
Young’s modulus mismatch, and the stress distribution in the oxide layer itself. Consideration of 
these factors leads one to the conclusion that delamination growth at this elevation into the oxide 
layer does not appear to be favorable. First, a weak bond coat/oxide layer interface would tend to 
deflect the delamination along the interface towards the trough region (where a compressive state 
of stress exists). This crack deflection is also driven by the lower value of the bond coat’s Young’s 
modulus relative to that of the oxide layer (for a high density oxide film). Further, the state of 
through-thickness stress in the horizontal cross sections through the considered crest elevation is 
not favorable to delamination growth. The through-thickness stress distributions are presented in 
Fig. 11 at several elevations through the crest. 

In the case of the thin oxide layer, the horizontal stress distributions are given at three elevations, 
namely: 0.5 fim below and above the bond coat/oxide layer interface; and 1.0 jim above the oxide 
layer/top coat interface, or at 23 = 639.5,640.5 and 642.0 /im (see Fig. 3(b)). As is observed, 
the initially tensile through-thickness stress directly below the bond coat crest (at 23 = 639.5 fim) 
becomes compressive immediately to either side of the crest within the oxide layer, and exhibits an 
oscillatory behavior with increasing 22 coordinate, going into small tension, and relatively greater 
compression. In the top coat region the through-thickness stress increases to a maximum value of 
63 MPA at 22 = 20 (im before it converges to the distribution observed previously in the absence 
of the oxide layer. The compressive region around the bond coat crest would thus tend to retard 
delamination growth at this elevation. At the higher elevation, 23 = 640.5 /zm, the magnitude 
of the tensile stress directly above the bond coat crest decreases relative to the lower elevation 
while the extent of compression to either side of the crest increases. However, the extent of the 
compressive region is smaller due to only one oscillation. A relatively high tensile stress ( 148 MPa) 
occurs a small horizontal distance from the crest at 22 = 18 /xra, which in fact is greater than the 
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stress directly below the bond coat crest at the lower elevation. However, it occurs within the oxide 
layer at the oxide layer/top coat interface and drops down to 61 MPa on the top coat side of the 
interface. This stress decays with increasing £2 coordinate, converging to the previous distribution 
around £2 = 25 /xm. The above suggests the possibility of delamination initiation in the top coat 
region to either side of the bond coat crest, with subsequent growth away from the crest arrested 
by the compressive zone above the trough region. At the highest elevation, £3 = 642.0 fxm, which 
is purely in the top coat region, the magnitude of the tensile through-thickness stress directly over 
the crest is further decreased slightly and the decrease of the stress with increasing £2 coordinate 
is substantially more gradual than at the lower elevations. In fact, the through-thickness stress 
never becomes compressive, reaching a minimum, then increasing somewhat and finally decaying 
to the distribution observed in the absence of the oxide layer. The above results indicate that the 
introduction of a thin oxide layer reduces the maximum tensile through-thickness stresses in the top 
coat in the vicinity of the crest where delamination initiation and growth are likely to occur relative 
to those when the oxide layer is absent. These results also suggest that delamination growth may 
involve coalescence of two delaminations initiated directly over, and to either side and below, the 
crest of the oxide layer/top coat interface. 

In the case of the thick oxide layer, horizontal through-thickness stress distributions are given 
at four elevations, with the first three located at the same elevations as in the thin oxide layer case 
and the remaining 1 fxm above the oxide layer/top coat interface, or at £3 = 639.5,640.5,642.0 and 
644.0 /xm (see Fig. 3(b)). The trends in the horizontal through-thickness stress distributions are 
similar to those observed in the case of the thin oxide layer, with a few notable exceptions. The 
large tensile through-thickness stress directly at the crest of the bond coat/oxide layer interface 
decreases with increasing elevation until it becomes compressive within, at and just above the oxide 
layer/top coat interface (as already seen in Fig. 10), in contrast with the preceding case where only 
tensile stress was seen. The magnitudes of the compressive through-thickness stress in the oxide 
layer immediately to either side of the crest are substantially greater relative to the thin oxide layer 
case, and they extend further to either side of the crest due to the greater thickness of the layer. 
The off-peak stress oscillations are also relatively smaller, with the tensile off-peak values smaller 
than the values directly at or above the crest (with the exception of the highest elevation). The 
largest off-peak tensile through-thickness stress (89 MPa) in the top coat occurs at the elevation 
£3 = 640.5 /im and the horizontal distance £2 = 22 /im, which is greater than the corresponding 
stress in the top coat in the thin oxide layer case at the same vertical elevation but at £2 = 18 fxm. 
The above observations suggest that, as in the preceding case, the potential for the initiation of 
delamination directly at the crest of the bond coat/oxide layer interface caused by the high tensile 
stress induced by the oxide layer is mitigated by the clamping action of the compressive stress to 
either side of the oxide layer. Further, the maximum tensile stress that occurs in the top coat 
above the oxide layer/top coat interface (at £3 = 650 /tm, see Fig. 10) is significantly reduced. 
This effectively moves the delamination initiation site to the off-peak location in the top coat region 
slightly above the crest of the bond coat/oxide layer interface. 

The increase in the oxide layer thickness also decreases the magnitude of the compressive 
through-thickness stress in the top coat region over the trough at horizontal elevations passing 
through the crest. This is more clearly shown in Fig. 12 where through-thickness stress distri- 
butions at the same elevation £3 = 640.5 fxm are given in the interval 20 < £2 < 50 fxm for 
the standard interfaces without and with thin and thick oxide layers. Increasing the oxide layer 
thickness is thus expected to facilitate delamination growth through the top coat region above the 
trough. 
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The last set of results demonstrates the effect of the oxide layer thickness on the transverse 
stress component (T 22 in the plane of the layered TBC system. This stress component arises due to 
the thermomechanical property mismatch between the individual regions regardless of the presence 
or absence of the rough interface and may potentially lead to the formation of vertical cracks. 
The magnitude of the transverse stress also depends on the geometry of the TBC system, and in 
particular on the substrate thickness. In the absence of a rough interface, the through-thickness 
(along the coordinate) distribution of the transverse stress exhibits abrupt jumps or discontinu- 
ities across the boundaries separating the substrate, the bond coat and the top coat because of the 
aforementioned thermomechanical property mismatch and also because this stress is not a traction 
component along this direction. In the absence of an oxide layer, the rough interface smooths out 
these abrupt jumps while the introduction of an oxide layer gives rise to large localized transverse 
stress discontinuities within, and in the immediate vicinity of, the oxide layer itself due to the 
layer’s large Young’s modulus. During the heating portion of the thermal cycle, the transverse 
stress component is initially negative in the substrate and positive in both the bond and the top 
coat such that the through-thickness integral of this stress component is zero due to the general- 
ized plane strain boundary conditions imposed on the TBC system in its plane. Due to the high 
relaxation rates, the transverse stress decreases to zero at the elevated temperature and at the end 
of the thermal cycle becomes positive in the substrate and negative in the bond coat and the top 
coat. During subsequent thermal cycles, the transverse stress decreases to zero in the three regions 
during the heatup portion and retraces its path during the cooldown portion as observed previously 
with the through-thickness stress. 

Figures 13 and 14 demonstrate the effect of the oxide layer thickness on the transverse stress 
distributions along the vertical cross sections (along the £3 coordinate) passing through the standard 
interface’s crest and trough during the initial heatup portion of the cycle at 30 sec and after 
cooldown at 630, 1260 and 1890 sec (end of first, second and third cycles). In the absence of an 
oxide layer, the transverse stress smoothly decreases as the rough interface’s crest is traversed from 
the bond coat into the top coat, becoming less tensile and less compressive in the top coat region 
away from the crest at 30 sec and after cooldown, respectively, Figs. 13(a) and 14(a). An initial 
increase in the magnitude of the transverse stress is observed as the bond coat/top coat interface 
is approached through the rough interface’s trough, followed by an initially rapid and then a more 
gradual decrease, Figs. 13(b) and 14(b). The presence of the oxide layer has a very localized effect 
on the transverse stress in both cross sections. During the initial heatup, Fig. 13, large tensile stress 
occurs within the oxide layer which increases with increasing oxide layer thickness. The region just 
outside the oxide layer is also affected in a manner that depends on the oxide layer thickness and 
the vertical cross section. In the case of the thin oxide layer at the crest location, Fig. 13(a), the 
region just below the bond coat /oxide layer interface sees a substantially higher transverse stress 
relative to the stress in this region in the absence of an oxide layer. The opposite is seen in the 
case of the thick oxide layer whose presence now lowers the transverse stress in the regions just 
outside the layer itself. At the trough location, Fig. 13(b), on the other, the initial increase in the 
transverse stress at the bond coat/top coat interface in the absence of the oxide layer is suppressed 
by the introduction of a thin oxide layer, while the transverse stress in the top coat at the oxide 
layer/top coat interface is enhanced. As the oxide layer thickness is increased, the transverse stress 
just below the bond coat/oxide layer interface is furthered lowered, with no additional increase 
above the oxide layer/top coat interface. The presence of the oxide layer has a greater effect in 
the trough where larger tensile transverse stresses are observed within the oxide layer relative to 
the oxide layer at the crest. It should be noted that the maximum values of the transverse stress 
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during the heatup portion of the thermal cycle do not occur at the same time as was the case with 
the through-thickness stress component. The critical times fall within the time interval between 30 
and 50 sec and depend on the location and the oxide layer thickness. The distributions presented 
at 30 sec, however, do provide a good representation of the effect of the oxide layer’s thickness on 
the magnitude of the transverse stress. Similar trends are observed at the end of the first as well 
as subsequent cycles, Fig. 14, except now the sign of the transverse stress is reversed. Further, 
the magnitudes of the compressive transverse stress in the oxide layer at the crest and trough are 
larger than the magnitudes of the tensile transverse stress at the corresponding locations during 
the initial heatup. 

The above results indicate that the magnitudes of the transverse stress within the oxide layer 
are substantially greater than the magnitudes of the through-thickness stress for the present TBC 
geometry. Should a vertical crack develop through the oxide layer during the initial heatup due 
to the large transverse tensile stress at the crest and trough locations, this crack may propagate 
into the top and bond coats due to the favorable Young’s modulus mismatch. Further, failure 
in compression at sufficiently large magnitudes of the transverse stress must not be ruled out as 
well. These observations point to the desirability of maintaining the oxide layer thickness below a 
certain level in order to eliminate the possibility of either a vertical crack formation or a compressive 
failure mode due to excessively large tensile or compressive transverse stress. The desirability of 
optimizing the oxide layer thickness is further supported by the examination of the remaining 
two stress components not discussed thus far. The out-of-plane normal stress (<r n ) also becomes 
progressively larger with increasing oxide layer thickness during the initial heatup, and may produce 
cracks within the oxide layer transverse to the xi axis. In fact, this stress component during 
the initial heatup is larger than the transverse stress. Its magnitude after cooldown, however, is 
relatively small. The inplane shear stress (0^3) also increases rapidly with increasing oxide layer 
thickness and attains maximum values at locations to either side of the crest. Its effect on potential 
delamination during cooldown is not clear, however, due to the clamping action of the through- 
thickness stress normal stress at these locations. 

4 SUMMARY AND DISCUSSION OF RESULTS 

As demonstrated by Koolloos and Houben [12], pre-oxidizing the bond coat to a thickness of 2 
\xvn enhanced the durability of a PS Z TBC system during furnace tests. The oxide layer consisted 
predominantly of a-A^Os. Failures occurred by spallation of the top coat upon cooldown. The 
results of our analysis shed light on the effect of the oxide layer thickness on the through-thickness 
(and also transverse) stress at the critical crest and trough locations of a rough top coat /bond coat 
interface. Such analysis can aid in optimizing the initial oxide film thickness for a specific TBC 
architecture and material system. 

Increasing the oxide layer thickness leads to dramatic increases in the magnitudes of the through- 
thickness stress at the bond coat/oxide layer interface and also within the oxide layer itself at the 
crest and trough locations after cooldown. The interfacial stress becomes more tensile and more 
compressive at the crest and trough locations, respectively. The increase in the tensile through- 
thickness stress, however, is offset by the compressive stress that develops immediately to either 
side of the crest, thereby providing a clamping action that should prevent delamination cracks from 
propagating across the oxide layer and into the top coat. The magnitude of this compressive stress 
increases with increasing oxide layer thickness. For a given oxide layer thickness, the through- 
thickness stress within the oxide layer at the crest decreases with increasing vertical distance from 
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the bond coat /oxide layer interface in the direction of the top coat, becoming less tensile at the 
oxide layer/top coat interface in the case of the thin oxide layer and compressive in the case of 
the thick oxide layer. This is accompanied by smaller magnitudes of the compressive stress in the 
locations immediately to either side of the crest at fixed elevations. Above the crest of the oxide 
layer/top coat interface, a small increase in the stress relative to that at the interface is observed in 
the case of the thin oxide layer. In the case of the thick oxide layer, a sign reversal occurs and the 
through-thickness stress becomes tensile again. Its magnitude, however, is substantially smaller 
than in the thin oxide layer case. 

The above results suggest that the thickness of the oxide layer has an influence on the manner in 
which delamination growth in the top coat will proceed, as is typically observed in plasma-sprayed 
TBCs. In the absence of the oxide layer, the maximum through-thickness stress (65 MPa) upon 
cooldown occurs at the crest of the bond coat/top coat interface and so the delamination is likely 
to initiate and proceed along this elevation. Its growth, however, should be somewhat mitigated 
by the change in the sign of the through-thickness stress with increasing horizontal distance from 
the crest. The introduction of a thin oxide layer moves the maximum stress (52 MPa) in the top 
coat at the crest location 2 /xm above the oxide layer/top coat interface, making it likely that 
delamination initiation and growth will proceed through this point and not through the crest of the 
oxide layer/top coat interface where the stress is 42 MPa (provided that the interface is sufficiently 
strong) . Since the thin oxide layer reduces the maximum through-thickness stress in the top coat at 
the crest location relative to a TBC without an oxide layer, delamination initiation at this location 
should require a greater number of thermal cycles. However, the thin oxide layer also increases the 
through-thickness stress to 61 MPa in the top coat to either side of the crest at an elevation below 
the oxide layer/top coat interface, potentially creating another delamination initiation site. The 
delamination growth at this and higher elevation should be mitigated by the compressive through- 
thickness stress away from the crest (in the top coat region above the trough) which has the same 
distribution as that in the absence of the oxide layer. 

The introduction of a thick oxide layer decreases the maximum through-thickness stress in the 
top coat at the crest location to 21 MPa which occurs 7 /xm above the oxide layer/top coat interface, 
further decreasing delamination initiation susceptibility at this elevation. However, this is offset by 
an increase in the through-thickness stress in the top coat to either side of the crest at the lower 
elevations, thereby potentially moving the delamination initiation site from the location directly 
above the crest to a location in the top coat to either side of the crest. Further, increasing the oxide 
layer thickness facilitates delamination growth between the rough interface’s crests by reducing the 
compressive through-thickness stress in the top coat above the trough region. 

The above discussion indicates that the oxide layer thickness should be optimized in a way that 
decreases the delamination initiation susceptibility in the top coat directly above, or immediately 
to either side of, the crest, while reducing the delamination growth potential between adjacent 
crests. Such optimization must take into account subsequent oxide film growth during exposure to 
elevated temperature. Another reason for optimizing the oxide layer thickness are the large through- 
thickness and transverse stresses that develop in the oxide layer during thermal cycling. These large 
stresses can lead to the degradation of the oxide layer itself, thereby producing additional failure 
initiation sites. The magnitude of transverse stress also depends on the TBC geometry, and in 
particular the substrate thickness which also must be taken into account in choosing on optimal 
oxide film thickness. 
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The results of our study can also be compared and contrasted with those presented by Freborg 
et al. [11]. These authors have modeled the oxidation of the bond coat layer, allowing the bond 
coat thickness to continuously decrease with exposure to elevated temperature with increasing oxide 
layer thickness, while in the present study the bond coat thickness and the oxide layer are fixed. 
However, the thickness of the oxide layer relative to the rough interface’s amplitude was relatively 
small. Further, the effect of the oxide layer growth on the resulting stress field should be small due 
to the high relaxation rates within the TBC constituents at the elevated temperature at which oxide 
growth was modeled. Any significant differences in the stress distributions between the two studies 
should thus be attributable to the model specimen geometry (cylindrical specimen vs flat plate), 
the waveform used to model the rough interface (single sine wave vs double sine wave), and the 
method of solution (f-e vs HOTFGM). The accuracy of the higher-order theory was demonstrated 
in previous investigations [16]. In the study of Freborg et al. [11], stress distributions are presented 
for a TBC without an oxide layer, and for a TBC with a growing oxide layer, at the end of 4 and 
51 thermal cycles. The thickness of the oxide layer after 4 and 51 cycles was approximately 0.3 
and 1.0 jum, respectively, which sire somewhat smaller relative to the waveform amplitude them our 
thin and thick oxide layers, respectively. 

First, both studies indicate comparable magnitudes and distributions of the through-thickness 
stress in the top coat region at the crest location in the absence of the oxide layer. Also, both 
studies indicate that the through-thickness stress after cooldown at the crest of the oxide layer/top 
coat interface decreases with increasing oxide layer thickness, and becomes compressive when the 
oxide layer becomes sufficiently large. However, the magnitude of the compressive stress after 51 
thermal cycles in the study by Freborg et al. [11] is substantially smaller than the magnitude 
of the compressive stress in our study for the thick oxide layer, which is likely attributable to 
the smaller oxide layer thickness relative to the waveform amplitude. Further, while Freborg et 
al. [11] predict that the stress at the trough of the oxide layer/top coat interface changes from 
compression to tension after 4 cycles with an additional, but relatively smaller, increase after 51 
cycles, our results indicate that the stress becomes less compressive (but not tensile) with increasing 
oxide layer thickness at this interface. This cannot be explained by the smaller relative oxide layer 
thickness in the Freborg et al. [11] study, and is most likely due to the differences in the waveforms 
and specimen curvature employed in the two studies. However, stress reversal in the trough region 
from compression to tension (smaller than in the Freborg et al. [11] study) is predicted to occur 
by our model in the top coat at a small distance away from the oxide layer/top coat interface in 
the case of the thick oxide layer (but not thin). This stress reversal also leads to lowering of the 
compressive through-thickness stress above the trough at the elevations passing through the crest, 
while in the study by Freeborg et al. [11], tensile stress was obtained at these elevations after 51 
(but not 4) cycles. 

Despite these quantitative differences, which are most likely caused by the differences in the 
relative oxide layer thickness, the waveforms used to model rough interfaces, and the specimen 
geometry, the qualitative trends predicted by the two models are similar. Our results thus also 
provide independent support to the hypothesis proposed by Freborg et al. [11] regarding the 
initiation and subsequent growth of delaminations during bond coat oxidation in plasma-sprayed 
TBCs with rough interfaces. In this hypothesis, the delamination initiation occurs at the crest 
of the bond coat/top coat interface (or at the off-peak location) and its growth between adjacent 
crests is facilitated by the growth of the oxide layer which changes the state of through-thickness 
stress from compression over the trough to tension. The results of our analysis suggest that the 
change of the through-thickness stress over the trough regions depends on the specimen geometry 
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and perhaps more importantly on the character of the rough interface. 

Finally, recent work by Brindley et al. [19] indicates that tailoring the thermal expansion 
coefficient of the bond coat in the region adjacent to the rough interface can double the TBC’s 
durability. This tailoring was accomplished through a fine dispersion of submicron alumina particles 
in the bond coat region adjacent to the interface. The higher-order theory employed in the present 
investigation is ideally suited to modeling and optimizing such advanced TBC systems through the 
already described subvolume discretization. Future work will address optimal grading of the rough 
bond coat/top coat interface using the concept proposed by Brindley et al. [19]. 

5 CONCLUSIONS 

The pronounced stress relaxation within the TBC constituents due to exposure to elevated 
temperature during uniform heatup-hold-cooldown thermal cycling changes the state of through- 
thickness stress at the crest and trough locations of rough bond coat/top coat interfaces in plasma- 
sprayed TBCs. At the crest, the initially compressive through-thickness stress during the heatup 
portion of the first thermal cycle relaxes to zero at the elevated temperature, thus becoming ten- 
sile upon cooldown. The opposite occurs at the trough. During subsequent cycles, the through- 
thickness stress at both locations decreases to zero during the heatup phase, and returns to the 
same state upon cooldown as after the first cycle. This sets up the crest location as the initiation 
site for delamination. 

Local variations in the details of the bond coat/top coat interface roughness tend to have 
a localized effect on the through-thickness stress. Increasing (decreasing) the trough amplitude 
of a rough interface relative to the crest amplitude increases (decreases) the magnitude of the 
through-thickness stress at this location, with substantially smaller changes at the crest location. 
These localized stress perturbations may be substantial, depending on the relative change of the 
trough (crest) amplitude with respect to the crest (trough) amplitude. The introduction of an 
aluminum oxide layer between the top and bond coats has a much more significant influence on the 
through-thickness stress at the crest and trough of the bond coat/oxide layer interface, producing 
substantially greater amplification in the magnitudes of the through-thickness stress after cooldown 
which increase with increasing oxide layer thickness. The delamination susceptibility at the crest 
of the bond coat/oxide layer interface, however, is mitigated by a high compressive stress in the 
oxide layer to either side of the crest, thereby preventing delamination growth into the oxide layer 
and the top coat. Further, the oxide layer’s effect on the through-thickness stress at the crest 
of the oxide layer/top coat interface is beneficial, lowering the through-thickness stress at this 
location with increasing layer thickness even to the point of putting this interface in compression 
(but not the region above it). Unfortunately, the tensile through-thickness stress in the top coat 
to either side of, and just below, the crest also increase with increasing oxide layer thickness while 
the magnitude of the compressive through-thickness stress in the top coat region above the trough 
of the oxide layer/top coat interface decreases, facilitating delamination initiation and growth at 
lower elevations. These opposing trends suggest that the oxide layer thickness should be optimized. 
This is further supported by the high magnitudes of the through-thickness and transverse stresses 
that develop in the oxide layer itself with increasing layer thickness. The oxide layer thickness 
optimization must therefore take into account the strength of the alumina layer (which depends on 
the deposition process), including interfacial strength, and subsequent growth. 

The major features of the residual stress evolution in TBCs with rough top coat/bond coat 
interfaces with an oxide film generated herein using the higher-order theory for functionally graded 
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materials are consistent with the results of finite-element analyses reported in the literature. In 
particular, the results of our investigation support a recently proposed hypothesis regarding delami- 
nation initiation and propagation in plasma-sprayed TBCs with evolving bond coat oxidation, while 
pointing to the importance of interfacial roughness details and specimen geometry in modeling this 
phenomenon. 
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7 APPENDIX 


In outlining the basic structure of the higher-order theory for functionally graded materials, 
distinction must be made between the global coordinates xi, X2, 23 shown in Fig. 2 , and the local 
inplane coordinates 2 , x^ ] used to designate position within each subcell (£7) of the generic cell 
(<?,r). The local and the global coordinate directions designated by the same subscript are parallel 
to each other. The subcell designation ( 0 j) identifies the subcell’s relative position within a given 
cell along the global 22, 23 coordinates. Since there are four subcells in each generic cell, with two 
subcells along each of the global 22, 23 coordinates, the indices 0 and 7 take on the values 1 or 2. 


7.1 Thermal Analysis 


Let the functionally graded parallelpiped shown in Fig. 2 be subjected to steady-state tem- 
perature or heat flux distributions on its bounding surfaces. Under these circumstances, the heat 
flux field in the material occupying the subcell (£7) of the ( q,r)th generic cell, in the region 


2 


OJ) 


— 


( 9 ) 


«( 7 ) 


X 


< must satisfy: 


dq 


OgH 

ai'-” 


= 0, 


03,7 = 1 , 2 ) 


(Al) 


The components qf 1 ^ 
according to: 


of the heat flux vector in this subcell are derived from the temperature field 



grr(P~t) 



(2 = 2, 3 ; no sum) 


(A 2 ) 


where k ^ are the coefficients of heat conductivity of the material in the subcell (£7), and no 
summation is implied by repeated Greek letters in the above and henceforth. 

The temperature distribution in the subcell ($7) of the ( q,r)th generic cell, measured with 
respect to a reference temperature T ref , is denoted by T^\ This temperature field is approximated 
by a second order expansion in the local coordinates as follows: 


r r(0l) _ J'ifh) I ~(P)rp(P~l) , =.h)rp(fil) , 1 ,n-(P)2 

1 — J (00) ' *2 ■ £ flO) +X 3 2 ( 01 ) " ) 'o(' 5X 2 


L (10) 


( 01 ) 


h iq)2 

n 0 


( 20 ) 




L (02) 


(A 3 ) 


where which is the temperature at the center of the subcell, and (m, n = 0, 1, or 2 with 

m + n ^ 2) are unknown coefficients which are determined from conditions that will be outlined 
subsequently. 

Given the five unknown quantities associated with each subcell (i.e., , ..., T^j) and four 

subcells within each generic cell, 20 N q N r unknown quantities must be determined for a composite 
with N q , and N r rows and columns of cells containing arbitrarily specified materials. These quanti- 
ties are determined by first satisfying the heat conduction equation, as well as the first and second 
moment of this equation in each subcell in a volumetric sense in view of the employed temperature 
field approximation. Subsequently, continuity of heat flux and temperature is imposed in an av- 
erage sense at the interfaces separating adjacent subcells, as well as neighboring cells. F ulfillm ent 
of these field equations and continuity conditions, together with the imposed thermal boundary 
conditions on the bounding surfaces of the composite, provides the necessary 20 N q N r equations for 
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the 20N q N r unknown coefficients in the temperature field expansion in each (/3q) subcell. 

The final form of this system of equations is symbolically represented below 


kT = t (A4) 

where the structural thermal conductivity matrix n contains information on the geometry and 
thermal conductivities of the individual subcells (/3q) in the N q N r cells spanning the and £3 
directions, the thermal coefficient vector T contains the unknown coefficients that describe the 
temperature field in each subcell, i.e., 


T = 


,(n) 

•11 


...,T 


( 22 ) 

NqN r 


where 

Igr ^ [ 2 ] 00) ? 2 T( 10 ) , T( 0 i ) , 2(20) > T(02) } ?r 7) 

and the thermal force vector t contains information on the boundary conditions. 


7.2 Mechanical Analysis 

Given the temperature field generated by the applied surface temperatures and/ or heat fluxes 
obtained in the preceding section, we proceed to determine the resulting displacement and stress 
fields. This is carried out for arbitrary mechanical loading, consistent with global equilibrium 
requirements, applied to the surfaces of the composite. 

The stress field in the subcell (/3q) of the ( q , r)th generic cell generated by the given temperature 
field must satisfy the equilibrium equations 




da 


dx 


*( 0 ) 


+ 


3 3 


dx 


(7) 


0, (j = 2,3) 


(A5) 


The components of the stress tensor, assuming that the material occupying the subcell (fyy) of the 
(q,r)th generic cell is isotropic, are related to the strain components through the familiar Hooke’s 
law: 


= 2m (^7) 
ij > 


Jfii) 


+ 


a(/37) 41 7) ^-2^ 7) 4 (/37) 


T(0 7) 

V 


(A6) 


where and (J® 1 - are the Lame’s constants of the material filling the given subcell (/5q), 

are the inelastic strain components, and aj^ 3 ^ is the thermal stress consisting of the products of 
the stiffness tensor components, thermal expansion coefficients and the temperature change. The 
components of the strain tensor in the individual subcells are obtained from the strain-displacement 
relations. 

The displacement field in the subcell (,/3q) of the ( q,r)th generic cell is approximated by a 
second-order expansion in the local coordinates x ^ , and £3^ 35 follows: 



= £i£n 
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where W^o}, which are the displacements at the center of the subcell, and the higher-order terms 
(* = 3) must be determined from conditions similar to those employed in the thermal 

problem. In this case, there are 40 unknown coefficients in a generic cell (g,r). The de- 

termination of these coefficients parallels that of the thermal problem. Here, the heat conduction 
equation is replaced by the two equilibrium equations, and the continuity of tractions and displace- 
ments at the various interfaces replaces the continuity of heat fluxes and temperature. Finally, the 
boundary conditions involve the appropriate mechanical quantities. 

Application of the above equations and conditions in a volumetric and average sense, respec- 
tively, produces a system of 40 N g N r algebraic equations in the unknown coefficients . The 

final form of this system of equations is symbolically represented by 

KU = f + g (A10) 

where the structural stiffness matrix K contains information on the geometry and thermomechan- 
ical properties of the individual subcells (^7) within the cells comprising the functionally graded 
composite, the displacement coefficient vector U contains the unknown coefficients that describe 
the displacement field in each subcell, i.e., 


U=[U< 1 1 1| ,...,U 


r (22) 
N q N r 


where 

= W m ,w i{1Q) ,W im ,W i{20) , w* 02)]g 7) (t = 2,3) 

and the mechanical force vector f contains information on the boundary conditions and the thermal 
loading effects generated by the applied temperature. In addition, the inelastic force vector g 
appearing on the right hand side of eqn (A10) contains the inelastic effects given in terms of the 
integrals of the inelastic strain distributions that are represented by the coefficients n j defined 
below 


R{ S?l,n) = ^’ A W /; /_; )Pn(& ] )dCf ^ (All) 


where 


A(mn) = ^\/(l + 2m) (l + 2n) 


and the Legandre polynomials P.(^) are functions of the non-dimensionalized variables £^’s, 
defined in the interval —1 ^ ^ ^ +1, which axe given in terms of the local subcell coordinates 
as follows: ^ /(h^ /2), and (3^ = 2). These integrals depend implicitly on the 

elements of the displacement coefficient vector U, requiring an incremental solution of eqn (A10) 
at each point along the loading path. 
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Generic cell (q, r) 


Figure 2.— A geometric model of a material functionally graded in the x 2 - x 3 plane 
illustrating the volume discretization employed in HOTFGM-2D. 
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Figure 3. — Volume discretization of the wavy interface region: (a) low, standard and high interfacial 
roughness; (b) standard interface roughness with thin and thick oxide layer thickness. 
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Figure 4. — Thermal loading history. 
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Figure 5. — Interfacial through -thickness stress continuity as a function of 
time across the crest and trough of rough interfaces without an oxide 
layer: (a) low interfacial roughness; (b) high interfacial roughness. 
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Figure 6.— Effect of interfacial roughness on the through-thickness 
stress evoiution with time at the bond coat/top coat interface: 

(a) crest location; (b) trough location. 
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*3 (fim) 

Figure 7. — Effect of interfacial roughness on the through-thickness stress 
distributions at the end of each thermal cycle (630, 1 260, and 1 890 seconds) 
along vertical cross section (in the X 3 direction) passing through the trough. 
Top coat/bond coat interface locations are given by X 3 = 626 pm (low 
roughness, dotted line), X 3 = 620 pm (standard roughness, dashed line), 

X 3 = 614 pm (high roughness, solid line). 
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Figure 8.— Effect of interfacial roughness on the through-thickness 
stress distributions at the end of each thermal cycle (630, 1 260, and 
1890 seconds) along horizontal cross section (in thex2 direction) 
passing through the crest of the bond coat/top coat interface. 
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Figure 1 0. — Effect of oxide layer thickness on the through-thickness stress 
distributions at the end of each thermal cycle (630, 1 260, and 1 890 seconds) 
along vertical cross section (in the X 3 direction) passing through the crest, 

X 3 = 640 pm (a) and trough, X 3 = 620 pm (b) of the standard interface. 
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Figure 1 1 . — Effect of oxide layer thickness on the through-thickness stress 
distributions at the end of each thermal cycle (630, 1 260, and 1890 seconds) 
along horizontal cross section (in the X 2 direction) passing below, through 
and above the bond coat/oxide layer/top coat regions at the crest of the 
standard interface: (a) thin oxide layer; (b) thick oxide layer. 
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Figure 1 2. — Comparison of the horizontal through-thickness stress 
distributions at the elevation X 3 = 640.5 p.m for the standard interface 
without and with thin and thick oxide layers. 
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Figure 1 3. — Effect of oxide layer thickness on transverse stress distributions during 
the initial heatup (30 seconds) along vertical cross section (in thex 3 direction) 
passing through the crest (a), and trough (b) of the standard interface. 
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Figure 1 4.— Effect of oxide layer thickness on the transverse stress distributions 
at the end of each thermal cycle (630, 1 260, and 1890 seconds) along vertical 
cross section (in the x 3 direction) passing through the crest (a) and trough (b) 
of the standard interface. 
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